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merically, including the determination of the parameters for univalence. The study 
involves the related mapping from an appropriate Euclidean rectangle to the circular- 
arc quadrilateral. Its Schwarzian derivative involves the Weierstrass p-function, and 
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Introduction 

Conformal mappings of the unit disk onto a symmetric circular arc quadrilateral 
(s.c.q.) with right angles at the vertices were investigated in [3llll[l5]. The space 
of conformal classes of such quadrilaterals is two-dimensional. One parameter 
may be taken to be the conformal module [1] of the quadrilateral, and the other 
as a real factor A appearing in the Schwarzian derivative of the conformal map- 
ping, which is a rational function whose general form is classical. The articles 
cited are principally concerned with the question of calculating these parameters 
numerically when given the specific geometry of the quadrilateral, and the range 
of values for which the Schwarzian derivative defines a univalent function. 

The aim of this article is to introduce an additional parameter in the problem 
by allowing a variable angle of avr radians at each vertex, where < a < 2. 
The numerical computations and the univalence criteria are generalized to this 
context and refined. 

An innovation here is to take into account the conformal mapping from a rectan- 
gle (rather than the unit disk) to the s.c.q. The advantage of this approach is that 
the Schwarzian derivative of the mapping takes the simple form (1/2— 2Q;^)p— 2//, 
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and the Weierstrass p-function has well known properties. The parameter /x G ffi 
plays the same role as the parameter A in [12]: in rough terms, either of these 
parameters controls the bending (i.e., curvatures) of the circular arcs bounding 
the s.c.q. We work out the explicit relationship between A and n by making use 
of the transformation (elliptic integral) of the unit disk onto a rectangle. Let 
± exp (±it) be the prevertices in the unit circle which map to the vertices of the 
rectangle; it is perhaps surprising that the expression for fi in terms of A depends 
on both t and a (see ffTB]) below). 

Our use of the p-function to clarify some aspects of the conformal mapping 
problem of s.c.q.s happily also sheds some light on the nature of the p-function 
itself. It is interesting, and apparently not well known, that the constants ei, 62, 
and 63 related to the p-function have simple expressions in terms of trigonometric 
ratios of t (see ffTOj) below). Further, we make some observations (in 12. 2p on the 
location of the zeros of p. 

It is well known that the conformal mapping of the unit disk onto a specified 
circular arc polygon can be expressed as a ratio of linearly independent solu- 
tions of a second order ordinary differential equation, with some parameters that 
have to be determined if the polygon has more than three boundary arcs. If the 
prevertices and the value of a are fixed, then the conformal mapping onto an 
s.c.q. has one unknown parameter. In [15] the second order differential equation 
is expressed as a Sturm-Liouville differential equation with parameter A, and 
the spectral parameter power series (SPPS) method is introduced to solve the 
differential equation. Here we use the SPPS method with parameter /i for com- 
puting the conformal mapping from a rectangle. A ratio of linearly independent 
solutions of the differential equation is in general (i.e., for any value of /i) con- 
formal locally, but need not be univalent. In Section H] below we show how the 
Sturm comparison theorem can be used to find bounds for n which are neces- 
sary for univalence of the conformal mapping. This makes it possible to find a 
good "seed value" for the SPPS method, and as explained and demonstrated in 
[T5] it is then possible to find numerically the exact bounds for fi for which the 
conformal mapping is univalent. 

As an application of the SPPS method (with a = 0) we compute the Fuchsian 
group for the universal covering map of a twice punctured disk. More precisly, 
for a unit disk punctured at points —a and a for any given value < a < 1 we 
are able to find the Mobius transformations of the Fuchsian group by calculating 
the value of a parameter s related to a. This is described in detail in Section 
below. As an application of the explicit calculation of conformal mappings, some 
plots are given of geodesies surrounding the two punctures for different values of 
a. 
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1 Description of problem 



In this article a symmetric circular quadrilateral (s.c.q.) is a domain D in the 
complex plane C bounded by a Jordan curve dD which is formed joining four 
vertices v, v, —v, —v by circular arcs with common internal angle an at each 
vertex (Figured]). By definition an s.c.q. is symmetric in the real and imaginary 
axes. 

1.1 Parameters of s.c.q.s 

First we need to describe the geometry of an s.c.q. D. Denote by pi, p2 the 
midpoints of the rightmost and the uppermost edges of D. Let 26i, 262 be the 
angles subtended by these edges with respect to their centers Oi, O2, and let Vi, 
r2 be the corresponding radii. Necessarily 



Mappings for a = 1/2 were studied in [3l HJ [15], and for a = in [2T] . 

Let pi, 6*1, and ri be given. To find v we need to know whether the right edge is 
convex. Let Ki, K2 denote the signed curvatures of the edges, that is. 
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Figure 1: Geometric parameters defining symmetric circular quadrilaterals. 
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so ri — l/\Ki\. If Oj = oo then Ki — while 9i is not defined. 

Let a also be given. Suppose first that Hi > 0. The tangent direction of the 
right edge at v leading towards pi is e**^^^~'^/^^ (necessarily pointing into the 
right half-plane), so the direction of the upper edge at v leading towards p2 is 
g«(ei-7r/2-a7r)^ It follows that the nature of bending of the upper edge is given 
by sgnfi;2 = sgn (sin(6'i — n/2 — an)) (although we cannot yet calculate ^2). If 
sgn K,2 < 0, then O2 lies on the ray from v in the direction of e*(^i~'^~"'^), whereas if 
sgn/t2 > 0, the direction is e**^^^""'^''. In either event this ray meets the imaginary 
axis in O2 = iilmv — Ret>tan(^^i — an)), and finally the values r2 — \v — O2I and 
then K,2 = (sgn/t2)/r2 are determined. 

Suppose on the other hand that Ki < 0. Then the tangent direction of the 
right edge at v leading towards pi is e~^^^^'^^^'^^ (necessarily pointing into the left 
half- plane), so the direction of the upper edge at v leading towards p2 is now 
g-i(ei+7r/2+a7r) ^j^gj^ sgu ^2 is the sigu of the imaginary part of this number, and 
O2 is calculated similarly. In summary, starting from 

1 

n -- 



1^1 



Oi = Pi - — , 

^ ^ (9^ + ^^g»(7r/2+(sgn«i)(ei-7r/2))^ 

one can calculate successively the parameters 

sgnK2 = sgn sin ( (sgn Ki) 6*1 - (1/2 + a)n), 

O2 — (Imt> — Rei'tan((sgn ft:i)6'i + an))i, 



K2 



Sgn K2 
\v-O2 



P2 = O2 + —i, 

«2 

02 = (sgn K2 ) arg — . 

V - U2 

If one knows only O2 and 62, there may be two possible upper boundary edges 
for D, which are complementary arcs of the circle centered at O2 and passing 
through V. Thus sgn K.2 is an essential datum. These formulas can also be applied 
when Ki — 0, with the convention that sgn = 0. 

We have shown that D is determined by pi, Hi, 9i, and a. Further, if pi and Hi 
are replaced by 1 and pi^i, the figure is rescaled. Thus for purposes of defining 
conformal equivalence classes we may normalize pi = 1 and conclude that there 
is a three-dimensional space of conformally distinct s.c.q.s. 

However, arbitrary values of the parameters given as above will not necessarily 
define an s.c.q. To begin with, ki and 9i must be such that the open right edge lies 
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entirely in the right half plane; for Ki > this requirement is pi — ri(l — cos ^^i) > 
0, while for ki < any value of 6i in (0, vr] is possible. Further, a must be chosen 
so that p2 does not lie in the lower half-plane. When ki > 0, the upper edge 
leaving v must not point simultaneously into the right half-plane and the lower 
half-plane, so if ki > 0, the largest avr can be is Stc /2 + 6i, in which case the 
upper edge of D goes horizontally from v to -|-oo and continues from — oo to —v. 
When Hi < 0, there are two cases. When a < 1/2, 6i may attain its maximum 
value of 77, in which case v coincides with the vertex v. When a > 1/2 the largest 
possible value of 9i is 37r/2 — avr, and again the upper edge of D goes horizontally 
through oo. We will loosely refer to these s.c.q.s as "of extremal type", when 
either t;=I;G]RorImt; = Imp2- 

1.2 Schwarzian derivative and accessory parameters 

Let D = {\z < 1} denote the unit disk. Let D be an s.c.q. and let /: D — )■ D 
be a conformal mapping such that /(O) = 0, /'(O) > 0. The preimages of the 
vertices of D are represented as ±e^** where < t < 7r/2. We will write when 
we wish to emphasize the topological quadrilateral whose boundary is with 
these four prevertices prescribed, in the context of conformal mappings from the 
disk to another topological quadrilateral. 

The Schwarzian derivative [18] Sf of f can be expressed as Sf{z) = Ra,t,x{z) 
where 

4(l-a^)V^o,tW-2AV'MW, (2) 
(cos2t)^'' - 2^2 -t- cos2i: 
{z^ -2{cos2t)z^ + iy ' 

2 sin 2t 
z^ - 2{cos2t)z^ + r 

(in the Appendix of [5T] the justification for general circle-arc polygons is worked 
out in detail; however due to a typographic error (1 — a)^ appears in place of 
(1 — a^).) Observe that the parameters a and A in Ra,t,\{z) separate nicely in 
the formula ([2]); this fact will be important in the sequel. 

A basic question is to investigate the set of combinations of accessory parameters 
{a,t,X) for which / is a univalent mapping. The first parameter is the easiest 
to relate to the geometry of the s.c.q. D, being simply l/vr times the internal 
angle at the vertices. We will see that relationship between t, A and the other 
geometric parameters (for instance k,i,6i) also involves a. It is convenient to 
do as much as possible with t alone, without considering a and A; this will be 
accomplished in the next section. 



^0,t(^) = 
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2 Mapping from a rectangle 



In this section we carry out a change of variable, writing f{z) = g{w) where 
g: R ^ D is a conformal mapping from a Euchdean rectangle R to an s.c.q. D. 
Taking the domain to be R instead of D naturally changes the behavior of the 
mapping at the vertices. In [20] such mappings from a rectangle were studied for 
a = but without consideration of their relationship to the corresponding maps 
from a disk. 

We will fix the prevertices of / choosing < t < 7r/2. For 2; G D define the 
complex elliptic integral 

(3) 

Note that the values Ui = E{1) and U2 = E{i) are respectively real and imagi- 
nary; it is well known that the image 

R = E(Dt) = {w: -ui < Rew < Wi, -\u2\ < Imw < \u2\} (4) 

is a rectangle with vertices at iwa, icJs, where = Ui + U2- We define 

g = foE-\ (5) 

the unique conformal mapping from R to D such that g{0) = 0, g'{0) > 0. By 
construction g takes vertices of R to vertices of D. 




2.1 The Schwarzian derivative and the p-function 

Our main task is to calculate the Schwarzian derivative of /, which by the Chain 
Rule [18] is equal to 

Sf = (Sg o E)E'^ + Se. (6) 

When / is continued analytically by reflection of z across edges of Dj (i.e., maxi- 
mal arcs of 9D — {prevertices}), the image value w = E{z) is reflected along edges 
of R. Even numbers of reflections return z to its original value and effect transla- 
tions of w. Since the Schwarzian derivative annihilates Mobius transformations, 
Se is a single valued function on the Riemann sphere and must therefore be a 
rational function. 

To describe Sg and Sf explicitly, let us introduce the auxiliary function 

if{z) = p{E{z)+ujs) (7) 
where p = pa;i,w2 is the Weierstrass p-function [23] with primitive periods 2uji, 

2U2. 
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As customary, write Cj = p^oJi), i = 1, 2, 3. These values are real and satisfy 

62 < 63 < 6i. (8) 

Let Q denote the quarter-disk Q = {z E 3: n < aigz < 3n/2}, considered 
as a topological quadrilateral with vertices (— 6**, — i, 0, — 1). Since p sends the 
subrectangle of R with vertices at (0, Ui, u^, IJJ2) to the lower half-plane, it follows 
from ([7]) that ip{Q) is also the lower half-plane. By applying the Reflection 
Principle one sees that must be a rational function of order 4 and so we can 
write it as 

ip{z) =ao + + + + —. 

Z — 6" 2; + 6 " 2; + 6" 2; — 6 " 

Clearly ip{-z) = p{E {-z) + ^3) = p{-E{z) - ^3) = p{E{z) + U3) = ip{z) and 
similarly v?(2) = ^p{z). From these symmetries we find that CI3 — — Cli, CI4 — — CI2, 
02 = — oT- Further, 

ao = lim (f{z) = p{E{oo) + cug) = p{uj3) = 63 

since E{oo) is a period. We have so far 

, -^^j^e (ai6**) + Re (016"**) 
= + ' .^-2(cos2t).^ + l • 

Since ip{0) = p^u^) = 63 we have Re (ai6~**) = 0. From this it follows that 
V?"(0) = 8Re(ai6**). On the other hand, differentiating the definition ([7]) and 
using the classical facts p'^Us) = 0, p"{uj3) = 2(63 — 61) (63 — 62) together with 
E'{0) = 1, E"{0) = we arrive at 

/ ^ (ei - 63) (63 - 62)^^ , . 

— 2(cos2t)2^ + 1 

Since (f is completely determined once t is specified, we are led to seek a formula 
for <f{z) which does not involve 61, 62, 63. By the symmetry of E with respect to 
reflections in the coordinate axes, we know that <f{—i) = 61, = 62, which 

reduces to 

61 = 63 + 4 sin^ t, 

62 = 63—4 cos^ t. 

The resulting relation 61 — 62 = 4, while perhaps surprising, is a consequence 
of the normalization we have implicitly applied to the periods of p. From the 
well-known fact 61 + 62 + 63 = we deduce that 

2 . ^ 

61 = 2--cos2t, (10) 

o 

2 

62 = —2 cos2t, 

3 

4 

63 = - cos 2t, 
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and ([9]) takes the definitive form 

cn(z) = - COs2t . (11) 

We observe tliat tlie lialf-period ui = K can be identified witfi tfie complete 
elliptic integral in the terminology of Jacobian elliptic functions [23], and that 
(es — 62) /{ei — 62) = cos^ t = fc^ is the modulus of the function sn(2;|fc^) which is 
related to p by 

p{z) = 62 + (ei - e2)/sn2(Vei - 62 z) (12) 
which in turn by <^ is equal to (-2 - (2/3) cos2t) + A/sn^2z. 

Next we observe that g extends by reflection along the edges of R and is also 
symmetric under reflection in the coordinate axes. From this it follows that 
when a; is a period of p, the value g{C, + u) is related to g{Q by a Mobius 
transformation. Therefore Sg is doubly-periodic with periods 2uji and 2uj2. A 
simple local calculation verifies that Sg has double poles at the vertices of i?, and 
since g converts the right angles there into angles of avr, we conclude that 

Sg{C)=^-^^p{C + 0O,)-2^ (13) 

for some constant /i G M. It is essential for us to see how /i depends on t and on 
the geometry of the s.c.q. D. From ([3]) we find that 

^ , , z cos 2t)z^ + (cos4t - 5)^^ + 2 cos 2t 
^^^^ ~ (24_2(cos2t)z2 + i)2 • 

From this and ([6]),(j7]), we obtain 

Sf = (^^^-2^)E'2 + 5^ (14) 

(2/3)(l-4«2)cos2t-2/i 
z'^ - 2{cos2t)z^ + 1 

2(cos 2t)z^ + (cos At - 2(1 - Aa^) sin^ 2t - 5)z^ + 2 cos 2t 
^ {z^-2{cos2t)z^ + iy 

after some manipulation. By comparison with ([2]) we discover that 

2(«'-l) , . X 

^i=— ^cos2t + 2Asin2t (15) 

3 

Thus for fixed a and t, the correspondence A •(-)■ /i is linear, and A, /i increase 
simultaneously inasmuch as0<t<7r/2. 

When one has generated a mapping / via parameters a, t, A, the corresponding 
parameters for g in (IT^ are obtained by keeping the same value of a, using /i 
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given by (fT5|) . and evaluating the integrals -E'(l), E{i) to determine the defining 
periods for p. On the other hand, when working with fll3p the data may be a, 
U2/UJ1 and /i. Then one may proceed as follows. 

We will say that a pair of half-periods Ui G M, W2 ^ is normalized when 
the corresponding values of the p-function satisfy Ci — 62 = 4. A given ratio 
r = U2/UJ1 G iM. is realized by a unique normalized pair (ci;i,c<;2)- To calculate 
it, one considers the p- function with (non- normalized) primitive periods 1, r 
evaluated at the half-periods and calculates r = y^p(T72y^^"p(r72)/4, and then 
takes Ui = r, U2 = rr. This may be done easily via Jacobian theta functions, 
using the formulas 

Pirn = Y^Mor+Mon 
Pir/2) = -^{Mor+M^n 

in the notation of [23], with the elliptic modulus q = e*'^^. 

We observe that with normalized periods, as r runs through the imaginary axis, 
the triple (61,62,63) = (2 — 63/2, —2 — 63/2, 63) traces the straight segment in 
R3 from (4/3, -8/3,4/3) to (8/3, -4/3, -4/3). 

Our point of view for the conformal mapping problem is that although the 
Schwarzian derivative ffT^ is a transcendental function, in some respects it is 
simpler to study than the rational function ([2]). 



2.2 The zeros of p expressed as elliptic integrals 

We digress to observe that this setting sheds some light on the location of the 
zeros of p. In [6j a formula for the zeros of p is given in terms of generalized 
hypergeometric functions. The starting point for the derivation of this formula is 
a formula in [7] which is an improper integral involving the well known modular 
forms Eq{t) and A(r). The latter formula may in turn be derived from the 
elliptic integral 

Mr) = ± r ^==£^= (16) 
Jo \/4:X^ - g2X - 93 

for the zeros Zo{t) of pr, where 4x^ — g2X — gs = 4(x — 6i)(x — 62)(a; — 63). 

We make the observation (not mentioned in [6] and [7]) that another formula for 
zq in terms of an elliptic integral is 



1 f^V^^ dv , , 

^o(r) = ± , / ; , (17) 

Vei-e2Jo ,/i _y2^/i_ (^a^^^2 
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This formula is an immediate consequence of the formula (fT2|) . We derive yet 
another formula for Zq as an elliptic integral: it follows from ([7]) that if 7r/4 < 
t < 7i/2 then the assertion p{zq) = can be expressed as ip{—iyo) = 0, where 

zo = E{-iyo) + UJI + UJ2, 

for some real number yo in the interval < < 1. Thus by flTTl) . 

4 4(sin2 2t)yi 

3'°'2^+yo' + 2(cos2t)yo2 + l-°' 

which implies 

2 -3 + cos^ 2t + (sin 2t)V9 - cos"^ 2t 

= 2^^;2^ • ^^^^ 

We therefore have 

'•^0 dx 



Zo = Ui+ U2 



I 



V + 2(cos2t)a;2 
2.o/(i+yg) 



UJI+UJ2--I , (19) 



1 -t;Vl - ismU)v^ 



By (HHD we find that 

Ay^ 4cos2t 1663 



{l + y^y cos2 2t + 2cos2t-3 (eg + 4)(3e3 - 4) ' 
Consequently, (fT9l) becomes 



Zo = Ui + U2 - - / , , = . (20) 



2 



By analytic continuation, the formulas f[T7j) and fl20|) are valid in general; i.e., 
p{zo{t) I r) = is valid for Zo{t) computed using f|T7|) or fl20l) with Imr > 0, 
since ei(r), e2(r), e3(r) and ;zo(t) are analytic functions of r. 



3 Calculation of the image quadrilateral as a 
function of the accessory parameter 

In this section we will carry further the method initiated in \i5\ by which geo- 
metric characteristics of the image s.c.q. are represented as power series in the 
accessory parameter A of ([2]) or /i of (fT3l) . This is a straightforward application 
of the recent Spectral Parameter Power Series (SPPS) method for solution of 
Sturm-Liouville differential equations pH], which we restate here briefly to make 
the presentation self-contained. 
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3.1 SPPS method 



We will work on the real interval < 2; < 1. One defines the sequence In of 
iterated integrals generated by an arbitrary pair of functions {qo,cii) recursively 
by setting Iq = 1 identically and for n > 1, 

In{z)= r In-liOqn-liOdC (21) 

where g„+2j = g„ for j = 1, 2, . . . 

Proposition 3.1 [13] Let ipQ and ipi be given, and suppose that yoo is a nonva- 
nishing solution of 

y'L + ^0 2/oo = Aoo^l Voo 

on the interval [0, 1]. Choose qo = l/y^, qi = ipivlo ^'^^ define X^'^\ X^"^ to be 
the sequences of iterated integrals generated by (go, qi) and by (gi, go) respectively. 
Then for each A G C the functions 

oo 

yi = yooY.{X-XoofX^''\ 

k=0 

oo 

k=0 

are linearly independent solutions of the equation 

y" + ^oy = X^iy (22) 

on [0, 1]. Further, the series for yi and y2 converge uniformly on [0, 1] for every 
A. 

Thus two linearly independent solutions of the second order ordinary differential 
equation (!22|) . which depend on the parameter A, are expressed as explicitly 
calculable analytic functions of A. 

3.2 p2 and ki 

We now concentrate on the geometric parameters p2 and Ki described in section 
[H where the s.c.q. D is normalized by pi = 1. Suppose that a,t,X are given. 
The formulation we give here for general a is analogous to the case a = 1 /2 as 
presented in [13]. Let / be the solution of Sf = Ra,t,\ in (recall ([2])) normalized 
by 

/(0) = 0, /'(0) = 1, r(0)=0, (23) 
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and suppose that f(JS)) is an s.c.q. similar to D, specifically /(D) = wiD where 
Wi > 0. Further, as is common, take yi to be a solution of (1221) in D with the 
normalization ?/i(0) = 1, ^^(0) = 0, where we use ipo^t, 4'i,t in place of ipo, i'l in 
Proposition 13.11 Then one may calculate the right edge midpoint Wi = /(I) via 
either of the expressions in 

Jo yi{ur yi{z) 

where y2 is the solution of fl22]) normalized by 2/2(0) = 0, ^2(0) — ^- Similarly, the 
upper midpoint W2 = f{i) can be calculated by integrating along the segment 
from to i, or more easily, by using i?o,7r/2-i,-A instead of Ra,t,x and integrating 
from to 1. We will refer to the corresponding solutions of this latter integration 
as y^, 1/2 • Now midpoint of the upper edge is seen to be 

P2 = ^ = (25) 

It was observed for example in [H [15] that the curvature of the right edge of 
/(D) = wiD is equal to yii^Y — 2yi{\)y[{l) . This can be verified by writing out 
the relationship \f{z) — Oi\\f{z*) — Oi\ = 1/ (curvature)^ which relates points 
z, z* symmetric with respect to the unit circle, and examining the limit as z — )■ 1, 
2;* — )■ 1. The result holds for zero curvature as well although the proof must be 
modified for that case. The curvature of the right edge of D is wi times the value 
for wiD, 

2/1(1)2/2(1) -2y2(l)yUl)- (26) 

Let a,t be fixed. Then how do ^2, ^^i depend upon A? Suppose we know a 
parameter value \^ such that the conformal mapping fa,t,\oo whose Schwarzian 
derivative is Ra,t,Xoo is univalent. (The question of finding such Aqo will be dis- 
cussed in section Hj). Then it follows that the solution 2/00 of (!22|) with A = Aoo and 
with the normalization 2/00 (0) = 1, 2/^(0) = is nonvanishing in D, and hence in 
[0, 1]. This 2/00 niay be used to define the power series 2/1 [A], 2/2 [A] of Proposition 
13. 1[ The notation [A] refers here to evaluation at z = 1 of the function whose 
parameter is A. Then fl23]) and fl26p give us 

MX) ^ «1 (27) 

yi[Aj2/2[AJ 

K,iX) = 2/i[A]2/2[A]- 22/2^2/1 [A] (28) 

For evaluation of P2(A) at a single A it one may prefer to evaluate the four power 
series on the right side of (1271) separately rather than to work out a single series. 
For solving equations of the form P2(A) = c, one may rewrite this equation as 
2/i[A]2/2[A] — c2/i [A]2/2[A] = 0, thus avoiding division of power series. 



12 



Now we consider the map g: R ^ D given by ([5]), defined in the rectangle R. 
Let Tji and 772 be the solutions of 



V'(C) + {{\ - a')p{( + 63) - /i) viO = (29) 
with normalizations (?7(0), -^'(0)) = (0,1) respectively, so 

because g'{0) = /'(O) = 1 (recall that E'{0) = 1 in section We have g{uji) = 
wi, g{uj2) = W2 so it is easy to express p2- As to the curvature, we note that for 
C G -R, the symmetry is 

\g(C)-OM2co,-C)-0,\ ^ 



curvature^ 



from which we find the curvature of the right edge of wiD to be 2?7i(a;i)?7^(a;i), 
which is a simpler expression than for the mapping / defined in the disk. Thus 

^1 = "^^^A^ViMii^i)) = 2r^2{ujMuj,). (30) 

Suppose that a and the normalized pair of half-periods uji,u:2 (or their ratio r) 
are given. Analogously to the case of mapping from D, if we have a parameter 
/X(X) for which g is univalent, we can form power series in /x 

K^{li) = 2r^2WM- (32) 

centered at /ioo, where [/i] , ''72 [/^] I'efer to the values obtained by integrating the 
differential equation vertically from to the endpoint u:2- 

From the above we are able to calculate ki, p2 numerically given the accessory 
parameters of the Schwarzian derivative, both for maps of D and of R. Further, 
if we assume that a, t or a, r are fixed and one of the values ki or p2 is specified, 
then we can approximate the parameter A or /i arbitrarily closely by truncating 
the corresponding power series to a polynomial of sufficiently high degree and 
calculating an appropriate real zero. Finally, when the full geometry of D is 
specified via a, ki, and p2, we can calculate (t. A) or (r, /i) by an iterative proce- 
dure in which a guessed value for the modulus t or r is corrected by solving for 
A or /i to achieve either one of ki or p2, and then adjusting to make the other 
value come out correct as well. There is no essential difference to the way this 
was done in [15] for a = 1/2, so we will not give further details here. 
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Figure 2: Extremal domains for a = 0.2, t = tt/6, with Amin = —0.479608 (left) and 
Amax = 1.30611 (right). Values near the vertices were plotted by extrapolation. 

4 Parameters for univalence 
4.1 Interval of univalence for fixed a,t 

Following the line of investigation of [3l IH [15], let us assume that a,t are given. 
We investigate the interval Amin < A < Amax of values of A for which the map 
/ = fa,t,x is univalent. A picture of the set of univalence pairs (t, A) corresponding 
to a = 1/2 was given in |15j . 

By the symmetry of s.c.q.s, a univalent / cannot have a pole in D. As an aid to 
the discussion we introduce the interval A^in < A < Aj^ax of values of A for which 
/ has no poles in ©. Thus 

Ajiiin — Amin < Amax ^ Amax" (33) 

Recalling the discussion of domains of extremal type in section 11.11 we first 
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suppose that a < 1/2. Then A = Amax produces an extremal domain when the 
vertices v = /(e**) and /(e~**) coincide on the positive real axis. When this 
happens, v = Wi + 2ri where as previously, Wi = /(I) and Ki = —1/ri < 0. 
From elementary geometry it is seen that the center of the upper edge of the 
s.c.q. D is O2 = W2 — ir2 and that arg(f — O2) = —air. From this it follows that 

\U!2\ — ^2 

smvro; = 

r2 

which says 

|w2|k2 — sinavr — 1 = 

since K2 = I/T2 > 0. By section [3] we know how to express the left hand side of 
this equation as a power series in 

A, which we may solve numerically for Amax- (The art of locating the zeroes of a 
polynomial with real coefficients is quite well developed and we will not discuss 
it here.) Similarly, the value A = Amin for which the upper and lower edges of 
D degenerate into full circles may be calculated by applying the same procedure 
using n/2 — t, —A in place of t,\. Once Amax (or Amin) is found, it is easy to 
calculate the conformal mapping. Examples are given in Figure [2] (a < 1/2) and 
Figure [3] (a > 1/2). 

The remaining case a > 1/2 is actually simpler. The upper and lower edges of 
D degenerate into horizontal rays meeting at 00 precisely when f{i) = 00, i.e., 
we have AJ^^^ = -^max- The degeneration of the right and left edges, /(I) = 00, 
corresponds to \^^^ = Amin- 



4.2 Estimates for A*^^^, A*^^^^.. 

There are certain advantages to working with the rectangle map g because the 
transcendental function Sg in many ways behaves more nicely than the rational 
function Sf. Thus we will investigate the values /imin^/^max which correspond 
-^min) -^max (flSjl . Siuce p IS real on the coordinate axes as well as on the 
boundary of the rectangle with vertices ibws, icJa, and since p' only vanishes 
at half-periods which are not periods, it follows from (fT3|l that the real function 
Sg is one-to-one on the horizontal segment from to uji and also on the vertical 
segment from to (X'2- This can be made more precise noting that p decreases 
from 63 to 62 on [u^, + Ui] while it increases from 63 to ei on [ci;3, + uj2]- In 
contrast, it appears to be much more difficult to determine bounds for Sf on the 
coordinate axes in 3. 

Lemma 4.1 Let + a{()ri{() = on a real interval [Co, Co + L], and let 

ri{Co) = 1; ''7'(Co) = 0. (i) Suppose a{C) < for all C G [Co, Co + -^]- Then rj never 
vanishes, (ii) Suppose a(C) > (vr/(2L))^ for all C- Then rj vanishes at least once 
in the interval [Co, Co + ■ 
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Proof. The proof is immediate from the Sturm comparison theorem [8]. For (i), 
note that the similarly normalized solution of f]"{C) + ■ ri{() = is = 1, which 
never vanishes. For (ii), use comparison with ^^"(C) + (^/(2-^^))^ • ^^(C) = 0, whose 
solution cos((7r/ {2L)){( — (q)) vanishes at (q + L. □ 

Proposition 4.2 Let the accessory parameters a and t (or equivalently t) be 
fixed, (i) If a < 1/2, then 



TT 



(ii) Ifoi> 1/2, then 



4 ' \2\uj2 



4 2 

,1 2\ / TT 



^4 ' \2\uj2 
Proof. From (fT3l) we are led to look a.t g = T]2/rii with r^i the solution of 

v"iC) + {i\ - a')p{C + CO,) - MO = (34) 

on the interval < < Ui, normalized as in Lemma 14. 1[ First suppose that 
a < 1/2. Then by ([HD the coefficient of rj in satisfies the bounds 

- a^)e2 - /i < - a^)p(C + ^3) - /i < - 0^)63 - /i, 

so rji must vanish at some point if — a^)e2 — /U > (7i^/(2wi))^, which tells us 
that fi < /i^jjj for such values of fi. Thus the first inequality of (i) is established. 
Further, we note that according to the Lemma, rji will never vanish if > 

(i-«>3. 

To examine the behavior in the imaginary direction it is helpful to introduce 
fi{u) = ri{() where ( = iu for < u < \u2\. The differential equation is 

ri"{u) + (/i - - a'^)p{iu + ^3)) ri{u) = 0, (35) 

so the normalized solution 171 must vanish if yU — (^ — a;^)ei > (7r/(2|u;2|))^, which 
tells us that fi > /i^ax such values of fi. Thus the final inequality of (i) is 
established. Further, we note that fji will never vanish if < (| — 0^)63. 

Part (i) has been proved. The proof of (ii) is similar; the main difference is that 
the inequalities derived from are reversed when 1/4 — is negative. □ 
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From two facts which we noted in the course of the proof, rj and fj do not both 
vanish for a common value of /i; i.e., g does not simultaneously have poles on 
both coordinate axes. The relationship of the value — 0^)63 to /iJnin, /^max; 
univalence appears to merit further investigation. 

From fl33|) we have immediately the following bounds on /imin, A'-max- 
Corollary 4.3 For a < 1/2, 



It follows from ([8]) that (1/4 - a^) 62 - {tc / {2uJi))'^ < /imin and /imax < (1/4 - 
a^)ei + (7r/(2|a;2|))^ for all a. The bounds can easily be converted to bounds for 
by means of ([T5]). 

Bounds for Amin, A^ax may also be derived in another way, by using Nehari's 
theorem [T7J which affirms that / holomorphic in D is not univalent when 



Thus a necessary condition for univalence is |5/(0)| < 6, but by f|T^ we have 
iS/(0) = 2(1 — a^)e3 — 2/i, so we arrive at the following result. 

Proposition 4.4 (1 — a^)e3 — 3 < yUmin and /imax < (1 — a^)e3 + 3. 

Numerically we can find values of a, t to make the bounds of either one of Propo- 
sitions 112] and 1131 better than the other (an example is shown in Figure H]). None 
of the above arguments can be expected to provide sufficient conditions for uni- 
valence, since our reasoning only samples a small subset of the domain of the 
mapping. 

A natural candidate for a value of A within the range of univalence is the average 
of the upper and lower bounds. Numerical experiments were made with the 
average of the bounds given by Propositions 14.21 and 14. 4[ as well as averaging the 
combined bounds (the maximum of the lower bounds and the minimum of the 
upper bounds). All of these appear to provide univalent mappings over a wide 
range of values of a and t, the best results coming from the Nehari bounds of 
Proposition 14.41 

Once one has succeeded in locating a value Aoo (or fiao) within the range of 
univalence, i.e., for which yi is a nonvanishing solution of the second order linear 




whereas for a > 1/2 




sup(l - \z\'^)\Sf{z)\ > 6. 



(36) 



17 



4 





0.5 1.0 1.5 2.0 



-2 



-4 



Figure 4: Comparison of upper and lower bounds obtained for /Xmim /"max for a = 0.25 
and 0.1 < r/i < 2.0, given by Proposition 14.21 (solid lines) and Proposition 14.41 (dashed 
lines). 

differential equation, tlie SPPS metliod enables calculation of Ki,Wi, etc., to 
considerable accuracy. (Actually, for the SPPS method to work it is not necessary 
to be within the domain of univalence, but only for the mapping to have no poles 
on the axis on which one is integrating.) Thus it is a simple matter to determine 
precise values of Amin, Amax numerically, as the nearest zeros of an appropriate 
power series (approximated by a polynomial) centered at Aqo- This was done in 
[15] for a = 1/2 and there is no essential difference for general a. 

In Figure [5] we show the A-domains of univalence for a = 0.0, 0.2, 0.4, 0.6, 0.8, 
and 1.0, together with the bounds described above. (The pictures for a > 1 turn 
out similar to a = 1) Following the presentation in |T5], we have applied the 
transformation arccot4A to make the domain finite. For a > 1/2 the bounds 
become difficult to distinguish from the contours of the domains. Domains for 
a > 1 are similar and are not shown. The median curves dividing the these 
domains of univalence in symmetric halves were devised as follows: The slopes 
of the boundary curves at the corners (0,0) and {tt/2,2tt) were estimated nu- 
merically based on the calculated values of Amin or Amax; write /3 for the slope at 
(0, 0). Then the median curve is given by the formula 



This appears nicer than a simple average of the upper and lower contours, and 
generalizes the dividing curve for a = 1/2 which is a straight segment as shown 
in [15]. One may use this median curve as starting values Aoo for the SPPS 
method. 
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arccot 4A arccot 4A 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 (1.2 0.4 0.6 O.S 1.0 1.2 1.4 



arccot 4A arccot 4A 




Figure 5: Values of arccot 4A for which the mapping fa,t,x (0 < t < 7r/2) is univalent 
in the unit disk. Bounds derived from AiminiA*max shown with a solid gray line; 
bounds derived from Nehari's theorem with a dashed line. 



5 Universal cover of a doubly-punctured disk 

As an application of the calculation of the accessory parameters of mappings of 
s.c.q.s, we consider the triply-connected plane domain = D — {—a, a} where 
the value < a < 1 is fixed. The universal covering map if: © — )■ has been 
studied in 01101122] • We can assume H{0) = and H'{0) > 0. The first question 
is to calculate or estimate H'{0). Another is to calculate the generators of the 
covering group of Ga corresponding to the covering H. We will use many basic 
facts about Fuchsian groups and covering spaces of Riemann surfaces, for which 
background information may be found in |16j . 

The mapping scheme is depicted in Figure |6l Write 

Ti{z) = — . 

z + I 
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Then Ti sends the upper half of D to its lower half, and when 

' = (tt5) • (^'> 

< t < 7r/2, we have Ti(e**) = a. It can thus be seen that Ti sends the upper 
half plane, punctured at the two points e** and — e~**, onto Ga- 

Now we consider the mapping problem of the disk to an s.c.q., taking a = 
and t as just defined. Let Amax be the parameter for which /(e**) meets /(e~**) 
as in section 14.11 where / = /o,t,Amax- Recall that this common vertex may 
be expressed as /(e**) = Wi — 2/ki (note that ki = — |ki| is negative) so we 
may divide by this value to obtain a mapping sending e** to 1 and sending 1 to 
l/{w,-2/K^). 

When we apply Ti to this rescaled s.c.q., the part in the upper half-plane is sent 
to a region in the lower half-plane which may be described as follows. It is the 
lower half of a domain Z)^ ^ D which is symmetric in the real axis and bounded 
by two arcs of dB) ending at ±e^*'^ and four arcs orthogonal to 83. The latter 
arcs join 1 to e^** and join —1 to — e^*'^, with tangencies at 1 and —1. By the 
above construction, we can deduce that 

e-'' = Ti( ^^^). (38) 

Finally, we define the conformal mapping H: Dg ^ Ga cls the composition 




Figure 6: Universal covering mapping of doubly-punctured disk. 



20 



Figure 7: Universal covering of doubly-punctured disk. 



From the above discussion, H sends the upper and lower arcs of (9D bounding 
the circular hexagon to the upper and lower semicircle of respectively, 
and sends the four remaining arcs into the segments [— 1, — a] and [a, 1]. Thus 
H may be extended by reflection to all of D, and together with one of its 
reflected images forms a fundamental domain for a Fuchsian group of the second 
kind [16j acting on ©, generated by two Mobius transformations Ai,A2 which 
preserve © and satisfy 



The quotient of D by this Fuchsian group is conformally equivalent to Ga, and 
the conformal mapping H (extended to all of D by reflection) is the covering 
map. 

The questions posed at the beginning of this section may now be answered as 
follows. Given a, we calculate t by fl37|) . and find Amax as in section I^?T1 Based 
on this parameter we determine Wi,ki as described in 13.21 This permits us to 
obtain s via (1551) . and then it is straightforward to find Ai and A2. Note that 
H{0) = 0, because 




e' 



±2, 
1, 



IS 



Alio. 



I, Ti(2) = 0, 



and since 



m)'(o) = 2, 



1 



(Tri)'(z) 



1 



f(0 



2 
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a 


t 


s/tt 


H'iO) 


tiAi o 


0.1 


1.37146 


0.32213 


0.00865 


11.027 


0.2 


1.17601 


0.26104 


0.58362 


19.171 


0.3 


0.98788 


0.21213 


2.04028 


31.388 


0.4 


0.80978 


0.17018 


3.73469 


51.328 


0.5 


0.64350 


0.13330 


5.31735 


86.582 


0.6 


0.48996 


0.10052 


6.61668 


155.776 


0.7 


0.34934 


0.07123 


7.57377 


314.856 


0.8 


0.22131 


0.04497 


8.19717 


797.411 


0.9 


0.10517 


0.02132 


8.53008 


3560.92 



Table 1: Numerical solution to universal cover of double punctured disk. (About 2.5 
seconds of CPU time to compute the entire table). 



the Chain Rule gives us 

In the notation of i^^, f'{i) = yi{i)~'^ which we have seen also can be expressed 
in terms of the power series produced by the SPPS method. Thus we may 
calculate H'{0). 

There is a unique simple geodesic in the hyperbolic metric of Ga which loops 
once around the points ±a. This geodesic is covered by the deck transformation 
Ai o of D, which fixes the points 

-1 + e'' ± V2 + 2e^'' 

on (9D. The trace of a hyperbolic Mobius transformation serves as a measure of 
the length of the invariant geodesic which it covers. These data can be calculated 
for any < a < 1 ; sample values are given in Table [H It is not difficult to see 
that as a — )■ 1, y4i o A^^ and the length of the geodesic tend to oo. However, 
lima_>.i H'{0) apparently is finite. It would be interesting to see if the values we 
have calculated could be obtained as limiting cases of the corresponding values 
for D minus two disks (pair of pants) as the radii of these disks tend to zero. 




Figure 8: Simple hyperbolic geodesic surrounding the two punctures of a disk. 
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It would also be interesting to determine the maximum of the crossing point of 
these geodesies with the imaginary axis (Figure E]). 

Conformal mapping of s.c.q.s can be useful for investigating many types of finitely 
generated Fuchsian groups. In general, such groups may define branched cover- 
ings rather than true universal coverings; a ramification of order n corresponds 
to a cycle of vertices [12] in the fundamental domain whose total angle is 27r/n. 
Thus it is of interest to map circle domains with diverse angles at the vertices. 
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